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FOREWORD 


The  flight  of  a  gliding  parachuta  in  a  uniform  wind  field  is  examined  using  the  methods 
of  Optimal  Control  Theory.  Thh  report  pres  nts  this  analysis  and  is  part  of  a  continuing 
effort  directed  toward  investigating  methods  which  will  improve  the  accuracy  and  dispersion 
characteristics  of  airdrop  systems. 

This  study  was  conducted  uruer  Department  of  the  Army  Project 
No.  1F162203AA33  Drop  Zone  Dispersion  Studies  by  Dr.  Allan  E.  Pearson,  Professor. 
C* i vision  of  Engineering,  Brown  University  while  he  was  assigned  to  the  Airdrop  Engineering 
Laboratory  as  c  visiting  scientist. 

Appreciation  is  extended  Mr.  Thomas  F.  Goodrich,  Mr.  Arthur  L.  Murphy,  Jr.,  and 
Dr.  Edward  W.  Ross  of  the  U.  S.  Army  Natick  Laboratories  for  their  help  Juring  this 
study. 
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I.  Introduction 

The  automatic  control  of  gliding  parachute  systems  has  been  considered  for  the 
recovery  of  sounding  rockets  and  in  various  military  airborne  operations.  Control  can 
be  effected  hy  a  servo-motor  pulling  on  the  shroud  lines  of  the  parachute  which  causes 
a  banked  turn  of  the  parachute  and  a  change  in  the  direction  of  flight.  Flight  tests 
of  an  automatic-homing  control  system  have  been  reported  [1]  which  detect  the  angular 
error  between  the  line  of  flight  arid  the  target  by  means  of  a  ground  based  transmitter 
and  a  pair  of  directional  antennas  located  on  board  the  parachute.  Referring  to  Fig.  1 
where  p  Is  the  position  vector  of  the  parachute  relative  to  the  fixed  target  in  the  horizontal 
plane,  v  is  the  horizontal  component  of  the  velocity  vector  of  the  parachute  relative  to 
air,  and  w  is  the  wind  velocity  vector,  the  automatic-homing  control  system  attempts 
to  align  v  with  -p,  i.e.  to  keep  the  angle  7  near  zero  at  all  times.  Goodrick  [2]  analyzed 
the  wind  effects  on  a  bang-bang  (full-on  or  full-cff)  version  of  this  control  system  and 
showed  that  Increasing  the  horizontal  system  velocity  relative  to  air  tends  to  decrease 
the  terminal  distance  error  from  the  target  in  the  presence  of  wind  distortion.  Murphy 
[3]  analyzed  the  performance  characteristics  of  the  automatic-homing  control  system  for 
the  ideal  situation  In  which  the  angle  7  in  Fig.  1  is  maintained  null  at  all  time.  Under 
the  condition  that  the  magnitude  of  the  wind  velocity  vector  w  be  less  than  the  magnitude 
of  the  system  horizontal  velocity  vector  v  relative  to  air,  i.e.  ||w||<||v||,  the  parachute 
possesses  a  wind  penetrating  capability  and  the  potential  of  reaching  the  target  under 
arbitrary  wind  directions.  The  assumption  is  made  that  the  initial  altitude  is  properly 
chosen  in  accordance  with  the  wind  angle  Lw  and  the  ratio  IMI/i|v||. 

In  this  report  the  control  of  a  gliding  parachute  is  viewed  in  the  context  of  optimal 
control  theory.  The  basic  philosophy  is  taken  that  at  some  time  tQ  intermediate  between 
the  launch  time  9  and  the  terminal  time  T,  0<  tg  <T,  an  estimate  &  of  the  wind 
vector  is  made  available  based  on  measurements  taken  In  the  interval 0<  t  <  tQ.  Under 
the  assumpt  on  that  the  wind  remains  constant  over  t0  <  t  <  T,  an  optimal  control  problem 
Is  formuiataO  which  minimizes  the  terminal  distance  error  from  the  target.  The  solution 
to  this  problem  Involves  integrating  a  system  of  four  nonlinear  differential  equations  with 
mixed  boundary  conditions,  i.e.,  a  two  point  boundary  value  problem.  An  approximate 
solution  to  this  problem  is  obtained  which  results  in  a  sub-optimal  feedback  control  law. 
Since  the  sub-optlmai  solution  is  a  closed-locp  control  law,  this  feedback  controller  can 
be  used  on-line  as  new  measurement!  of  the  wind  vector  w  become  available.  The  filtering 
problem  of  estimating  the  wind  based  on  available  measurements  is  not  considered  in 
tnls  report. 


The  equations  of  motion  and  optimal  control  system  formulations  are  presented  in 
Section  II.  The  sub-optimal  control  solution  is  presented  in  Section  III.  Various 
extensions  of  this  work  are  indicated  in  Section  IV. 


II.  Optimal  Control  Problem  Formulations 

t 

Assuming  a  constant  wind  w  in  the  horizontp'  '•line,  a  constant  rate  of  descent  v  , 
and  an  initial  altitude  h0  at  launch  time,  the  equation  of  motion  governing  the  parachute 
can  be  represented  in  the  horizontal  plane  according  to 


p  ■  v  +  w,  0<  t  <  T  c  ^q/v  (1) 

In  this  report  the  velocity  vector  v(t)  will  be  assumed  to  be  constant  in  magnitude 
IMt)||  -T|v,(t)|J  +  |Vj(t)|*]  1/4  -  a  ■  constant. 

Then  v(t)  can  be  represented  by 
v,  (t)  ■  a  cos  co(t) 

(2) 

v3(t)  »  a  sin  fo(t) 

where  the  angle  <o  «  is  related  to  the  bank  angle  0  of  the  parachute  via  the  well  known 
relation 


co  a  tan  0.  (3) 

llvll 

Since  the  bank  angle  0  can  be  directly  manipulated  by  changes  in  the  servo-motor 
connacting  the  shroud  lines,  it  is  sufficient  to  summarize  (1)  -  13)  by 

p,  *  a  cos  to  +  w, 

p,  *  a  sin  co  +  wj  (4) 

co  *  u  to<t<T 

when  u  is  regarded  as  the  control  variable. 
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I 
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Let  p(to),  «(to)  and  w  be  given  at  some  Initial  time  ^  in  the  interval  0<  tQ  <  T. 
A  performance  index  which  reflects  a  number  of  desirable  features  for  this  problem  is 

P  -  H||p(T)||J  +  -SL|  w(T)  -  Z.w  _  ,13  +  _SL_  /]u(t)|2dt  (5) 

2  2(T-t0)  t0 

where  qt  and  qa  are  non-negative  weighting  parameters.  The  first  term  in  (5),  1/2||p(T)||J, 
reflects  tha  desirability  of  minimizing  the  Euclidean  distance  from  the  target  at  the  terminal 
tima  T.  The  second  term  with  weighting  parameter  q,  reflects  the  desirability  of  having 
the  parachute  point  upwind  at  tha  terminal  time  in  order  to  reduce  the  total  horizontal 
velocity  at  touch-down.  The  third  term  with  weighting  parameter  q2  reflects  the  cost 
of  control  effort  in  terms  of  the  “a  wage  power"  expended  over  the  interval  ^  <  t 
<  T. 


Before  proceeding  with  the  necessary  conditions  for  a  minimum  of  (5),  some  further 
comments  regarding  the  role  of  the  weighting  parameters  is  in  order.  In  the  absence 
of  any  other  constraints  on  tha  control  variable  u,  it  is  necessary  to  choose  qa  >  0  in 
order  to  be  assured  of  a  meaningful  solution.  Physically,  an  optimal  control  solution 
with  qa  small  will  laad  to  relatively  smaller  values  of  the  te'/ninal  distance  error  ||p(T)|| 
but  will  generally  require  relatively  larger  bank  angles  <p  of  the  parachute,  in  contrast 
with  larger  values  for  qa .  In  the  Interest  of  making  ||p(T)||  as  small  as  possible,  a  meaningful 
solution  doe7  result  for  qa  -  0  provided  a  constraint  on  u  is  assumed  of  the  form  |u(t)KM, 
i.e.  the  bank  angle  <t>  is  affectively  limited  by  v  ie  upper  bound.  However,  the  analytical 
solution  to  this  problem  is  less  tractable  than  t.  •  unconstrained  problem  with  qa  >  0. 
As  will  be  Indicated  in  Section  IIB,  the  solution  vi  the  constrained  optimization  problem 
is  bang-bang  followed  by  a  possible  null  interval,  i.e.  u(t)  *  ±  M  or  0,  which  implies 
instantaneous  reversal  of  bank  angles.  Since  the  dynamics  of  the  servo-motor  and  bank 
angle  rates  are  not  reflected  in  the  equations  of  motion  (4)  the  value  of  this  problem 
(qa  -  01  is  probably  of  lass  practical  importance  although  it  is  of  theoretical  interest. 

A  further  comment  regarding  the  performance  index  (5)  refers  to  the  second  term 
with  weighting  parameter  q( .  There  might  be  some  motivation  for  choosing  q,  *  o  and 
adding  the  terminal  constraint  that  w(T)  *•  £w  +  n.  However,  as  in  liie  case  with  potential 
terminal  constralnst  on  p(T),  it  will  not  generally  be  possible  to  satisfy  such  constraints 
in  the  fixed  Interval  (t^  T)  tor  any  but  a  small  set  of  initial  conditions.  In  a  practical 
implementation  of  an  optimal  control  law,  q,  would  probably  be  chosen  as  zero  with 
a  programmed  command  that  with  A  T  seconds  to  go  die  parachute  execute  a  rapid  turn 
into  the  wind  regardless  of  the  distance  from  the  target. 


A.  Transformation  to  Normallzad  Variablas 

Let  e  time-varying  transformation  of  the  origin  be  made  according  to 
y  -  p  +  (T-t)  w  ,  tp  <  t  <  T. 

Then  minimizing  the  terminal  distance  ||p(T)||  is  equivalent  to  minimizing  ||y(T)||.  In 
addition,  let  the  independent  variable  be  transformed  via 


and  define  new  dependent  variables  via 


Yi 


a(T-to) 


,  x2 


Ya 


a(T-to) 


x3  *  to 


(8) 


In  terms  of  these  variables,  the  equations  of  motion  (4)  become 


x;  *  cos  x3 


Xj  •  sin  x3 


(9) 


x3  -  (T-to)u  (Kr<1 


where  prime  denotes  differentiation  with  respect  to  r- 


The  performance  index  (5)  becomes 
aJ(T-to)J 


P  - 


fc,(1)P  +  lxa(1)lJ  +  Q,|x3(1)-/w 

1  1  1 

•  i  A  8  L.  12  V 


-/W  -  rt3 


where 


+  Qj  /-Ju(r)pdr| 


(10) 


n  ■  ft  ■ 

°‘  aJ(T-to)J  ,  osrr-to):i 


(11) 
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Hence  minimizing  (5)  or  (tO)  itequiveient  to  minimizing 


Combining  (9)  with  (14)  -  (16),  and  atfining 

X4  -  X3 ,  (17) 

it  foiiows  that  the  necessary  conditions  for  the  unconstrained  optimai  control  probiem 
involve  solving  the  two  point  boundary  value  problem: 


x,’  “  cos  x3 


Xj  ■  sin  xs 

xj  *  -  X4 

Qa 

x4  *  x,(1)  sin  x3  -  x2(1)  cos  xj 

subject  to  the  given  initial  conditions 

x,(0)  “  - 1 -  (Pi  (to)  +  (T-to)w,) 

a(T- to) 

x?  (0)  - J —  [pj  (to)  +  (T-to)w,] 

a(T-t0» 

xj(0)  «  cj  (t^) 
and  the  terminal  condition 


(18) 


(19) 


x4(1)  -  Q,  [x3(1)  -  P  -  *]  .  (20) 

By  u  reduction  of  (18)  to  three  differential  equations  involving  dx,/dx3,  dx2/dx3 
end  dx4/dxSf  it  is  possible  to  derive  three  algebraic  equations  in  the  three  unknowns 
x,<1),  xa(1)  and  x^O)  which  depend  on  the  given  initial  conditions  in  (19)  However, 
these  equations  involve  Eiiiptic  integrals  of  the  first  and  second  kind  with  arguments 
depending  on  the  mixed  boundary  data  (see  Appendix).  Ht'ce  an  exact  analytical  solution 
to  (16)  -  (20)  via  this  approach  does  not  appear  promising,  although  some  approximations 
to  these  agnations  might  be  fruitful.  This  has  not  been  explored  in  any  detaii. 
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in  the  case  of  the  constrained  optimal  control  problem  in  which  Q2  «  0  and  u  is 
constrained  via  |u(tM  <  M,  the  Hamiltonian  becomes 

H(X,x,u)  *  X,  cos  x3  +  Xj  sin  x3  +  (T  -  t^  X3u 

where  (X,  ,Xa,X3)  satisfy  the  same  differential  equations  as  in  (14).  However,  in  this  case 
the  minimum  principle  leads  to 


u# 


The  appropriate  value  for  u*  over  a  singular  interval  during  which  X3  *  0  can  be  resolved 
by  the  condition  that  X3  «  0  ■»  X3  *•  0  x3  *  constant  ■*  x3  ■  0  -*  u*  "0.  Moreover, 
if  a  singular  control  Interval  Is  achieved  on  some  time  interval  contained  within  the  given 
interval  to  <  t  <  T,  the  control  remains  singular,  i.e.  null,  throughout  the  remainder 
of  the  Interval.  This  fellows  because  u*  *  0  •*  x3  ■  0  which  means  x«  remain?  at  zero 
for  the  remainder  of  the  in* erval.  Hence  the  optimal  control  is  bang-bang-c rf  for  the 
constrained  problem  with  no  intermediate  singular  intervals.  It  is  difficult  ro  say  how 
man*/  switches  will  occur  during  the  interval  tQ  <  t  <  T,  although  presumably  more 
switches  will  occur  when  the  parachute  is  close  to  the  target  with  a  large  amount  of 
time  to  go  until  touch  down. 

In  summary,  the  solution  to  the  constrained  optimal  control  problem  involves  solving 
the  two  point  boundary  value  problem: 

x,’  ■  cos  x3 

xi  ■  sin  x3 

x3’  -  -(T  -  to)  M  Sgn  (  4) 

xi  *  x,(1)  sin  x3  -  Xj(1)  cos  x3 

subject  to  the  initial  and  terminal  data  in  (19)  and  (20),  where  tne  "Sgn"  function  is 
defined  by 


C-M  if  X3  >  0 
J+M  If  X3  <0 

'  nn  conclusion  If  X3 


L 


-  0 
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Sgn(x) 


+1  If  x  >  0 
0  if  x  -  0 
-1  if  x  <  0 


The  open-loop  optimal  control  solution  u*(rp(t0),c>>lt0)),  0  <  r  <  1,  regarded  as 
a  function,  of  the  initial  conditions  p(t0)  and  w(t0),  becomes  a  closed-'oop  or  feedback, 
control  law  for  either  of  die  above  problems  simply  by  choosing  r  *  0  and  using  the 
continuous  variables  j.(t)  and  w(t)  in  place  of  p(tg)  and  w(t0).  in  this  way  u*»0,p(t),w(t)) 
is  the  optimal  control  for  every  t  based  on  continuous  measurements  in  p(t)  and  w(t). 

II!.  A  Sub-Optimal  Feedback  Control  Lew 

In  the  interest  of  obtaining  a  feedback,  i.e.  ciosed-loop,  contro'  law,  an  approximate 
solution  will  be  developed  in  this  section  to  the  necessary  conditions,  (18)  -  (20),  f  >r 
the  unconstrained  optimal  control  problem.  Let  a  parameter  €  be  defined  by 

g-  (T-to)*  -  aa(T-t°)4  (21) 

Qs  ds 

and  consider  (18)  -  (20)  in  vector  notation: 


x*  ■  f(x,r,€) 


cos  x, 


sin  xj 


Xj  (1 )  sin  xj  -  x,(1)  cos  Xj 


A  x(0)  +  B  x(1)  -  f 


where  the  matrices  A,  B  and  f  are  given  by 


10  0  0 
0  10  0 
0  0  10 
0  0  0  0 


0  0  0  J 

0  0  0  0 
0  0  0  0 
0  0  -Q,  1 


A 


and  the  constants  {afi.yfi)  are  defined  as 

a  «  x,  (0)  ■  — - - (Pi  (to)  +  (T-to)W|  ] 

a(T-to) 

0  -  x2(0)  -  - - —  (Pj(to)  +  (T-to)wj) 

a(T-to)  (25) 

7  *  Xj(0)  "  w(to) 

C  «  &  +.  ir 


In  view  of  the  analytic  properties  of  the  function  f(x,r,6),  a  Taylor  series  expansion  of 
the  solution  x(t,€)  to  (22)  and  (23)  ebout  6-0  will  exist  for  sufficiently  smai:  6.  Thus 


x(r,€) 


I  3nxin=)  I  €n 
n-0  aen  |  n! 

6-0 


(26) 


and  an  approximate  solution  to  the  problem  is  reoresented  by  retaining  a  finite  number 
of  terms  In  the  series  (28).  The  corresponding  epproximation  to  the  optimal  control 


u"  given  by  (15)  and  (17)  is  obtained  from 


T-to  n-0  3€n  €«0n! 


(27) 


The  sub-optimal  feedback  control  law  follows  by  choosing  r  >  0  in  (27)  and  letting  the 
initial  time  t<,  in  (25)  be  replaced  by  running  time  t. 

Before  proceeding  with  the  calculations  of  the  coefficients  in  (26),  it  should  be  pointed 
out  that  €  -  0  corresponds  ito  no  control  effort  (u  -  0)  expondod  on  the  interval 
tj,  <  t  <  T.  This  Is  evident  since  6-0  implies  Qa  -  <»  for  which  u  -  0  is  the  optimal 
solution  in  minimizing  (12).  In  other  words,  thv  nominal  solution  in  this  approximation 
is  the  uncontrolled  descent  of  the  parachute. 

In  order  to  simplify  the  notation,  let  the  components  of  the  vector  coefficients  in 
(26)  be  represented  by  (an,  bn,  c^dn)  according  to 


9 


anXj  tr,€) 
dcn 


€<«0 


*nM 

cn(r) 


anx>(r,e) 

ae  n  1 

€•0 


bn(r) 


anx3(T,e)  | 
aen  > 
€■0 


The  zeroed  order  <n  *  01  term  ».  (26)  is  then  obtained  by  solving 
a^  -  cos  Cq 

bo  -  sin  Cq 

-0 

do  *  aD(1)  sin  ®o  “  b0(1)  co*  °o 

subject  to  the  boundary  conditions 

ao(0)  -  a  ,  bo(0)  -  0  ,  Cq(0)  -  y 

d0(^)  -  Q,  [co(D  -  6) 

The  solution  to  (28)  and  (28)  is  readily  obtained  as 
ao(r)  "  a  +  tC 

b0(r)  -  0  +  tS 

Cp(r)  -  T 

d0(r)  »  (1 -r)D  +  Q,(7-«) 

where  C  -  cos  S  *  sin  y  and 
D  *  0C  —  aS 

t 

v 

have  been  used  for  shorthand  notation. 


(28) 
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The  first  order  (n  -  1 )  term  In  (26)  Is  obtained  by  differentiating  both  sides  of 
(22)  with  respect  to  €  and  then  putting  €  -  0.  This  yields  the  system  of  equations 
for  (a,,  b1#  c,,  d,): 


a,'  -  -Sc, 

bi  *  Cc,  (32) 

c,'  -  -d0 

dj  -  a,(1)S  -  bt(1)C  +  [«o<1)C  +  b0(1)S]c, 


The  boundary  conditions  for  (32),  as  for  all  subsequent  terms  in  the  series  (26),  are  given 
by 


36n  I 

eio 


a"x(i,6) 

d€n 


-  0 


6-0 


n>1. 


(33) 


Yhls  follows  from  (23!  and  the  fact  that  (A.BJ)  to  not  depend  on  6.  Thus  (33)  yields 


■n(0)  -  0.  bn(0)  -  0,  c„(0)  -  0 

(34) 

dn(1)  -  Q,cn(1K  n  >  1 

The  integration  of  (32)  subject  to  (33)  with  ag,  bo  and  do  given  by  (30)  yields 

a, (r)  -  -  S[-£  r*  +  (Q,  (y-3)  -  D)  1*  ] 

0  2 

b, (r)  -  C[fi  H  ♦  (Q,(7-«)-0)L*I 

6  2 

c,  (r)  -  (Q,  (y-«)  -  0)  r  ♦  ^  r*  (35) 

d,  (r)  -  d,  (0)  -  [HQ,  (7-«)  -  1/3  D) 


■Ml  +aC  ♦  ?S)[(Qj(7— 8)-0)  ~  ♦  D  ~) 

2  o 


ii 


% 


where  the  initial  condition  on  d,  is  given  by 

d,(0)  -  1/3(aC+j3S)  D  +  Q,  [(7-  6 - 1/2 D)Qt  -  1/zkx  +I80M . 

The  firv*  two  terms  of  the  sub-oftimal  feedback  controi  lew  (27)  with  t  -  0  end 
N  *  1  have  bet^  obtained  up  to  this  point  with 

M(0,e )  - - [do(0)‘  +d,  (0)€  1 

T-to 

■ - [1  +  — (ac+^S)]  D  +  (36) 

T-to  3 

[7*-C  +  e (7  —  5 - 1/2 D )Qj  (aC+/3S)}  Qj 

'-<o  2 

The  higher  order  terms  in  the  series  (26)  become  i  ncreasingiy  more  tedious  to  obtain. 
For  r  •  2,  fie  differential  equations  to  be  soived  'after  some  simplification)  are 

a, '  -  — C(ci  ]  *  — •  Sc* 

b, '  -  — S[Cl  1 2  +  Co, 

ci  -  -2d, 

d,  -  Saj(1)  -Cbj(1)  +  Die,!2  h-  [1  +  oC  +  /SS]ca 

where  c,  (r),  Cj(r),  d,  (r)  are  given  by  (35)  and  the  boundary  conditions  by  (34).  These 
equations  were  integrated  under  the  special  case  Q,  -  0  with  the  result  that  the  desired 
initial  condition  on  da  is 

Q,  -  0:  da(0)  -  D  [-1  -  ^D3  +  1  («C+0S) 

1  o  ID  40 

+  £  (aC+flS)*  ] 

15 

The  interesting  feature  about  the  case  Q,  ■  0  is  the  observation  that  the  sub-optimal 
control  law  contains  0  as  a  factor  -  at  (east  as  far  «  the  first  three  terms  in  (27)  have 
been  oomputed.  Thus  it  would  appear  as  though  the  quantity  0  in  (31), 

D  •  xa(0)  cos  Xj(0)  -  x,(0)  sin  x,(0)  (37) 
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plays  a  major  role  in  the  optimal  feedback  control  law. 


it  is  instructive  to  refer  the  above  obtained  sub-optimai  control  law  back  to  the 
original  coordinate  system.  Assuming  Qt  -  0  ir.  (36)  and  referring  to  (37),  (25),  (21), 
(9)  and  (7),  this  results  in 

u  -  -U= y'  Ev  M  +  -H=*L3  /v  +  . . . )  (38) 

<fc  .  3qj 

in  the  above,  y  ■  y(t)  is  the  moving-origin  coordinate  system  defined  in  (6),  i.e. 
y  -  p  +  (T  - 1)  w, 

% 

and  E  is  the  2  x  2  matrix 

■t:;) 

which  effectively  generates  the  D  term  in  (37)  ,  i.e. 


y'Ev  -  YiVj  —  YjVj- 


IV.  Concluding  Remarks 

Given  an  estimate  of  the  w'nd  vector  w  at  some  initial  time  to  intermediate  between 
launch  and  touch  down,  and  given  the  initial  position  and  velocity  vectors  (relative  to 
air)  of  the  parachute,  th>  optimal  control  of  the  parachute  has  been  considered  which 
strives  to  minimize  the  terminal  distance  of  the  parachute  from  the  target  at  the  time 
of  to  ich  down.  A  sub-optimal  feedback  control  Jaw  has  been  obtained  which  can  be 
used  in  a  closed-loop  manner  based  on  continuous  measurements  of  the  position  and 
velocity  vacrur*  of  the  parachute.  The  first  order  approximation  to  this  control  iaw  is 
strikingly  simple  and  takes  the  form 

(y,v2  -  yav,j  (30) 

3a 
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■f. 


where  Vi  ■  Pi  +  (T  -  t)wi  and  y2  ■  p2  +  (T-  t)w-  are  the  horizontal  component* 
of  the  position  vector  p  -  p(t)  of  the  parachute  at  time  t  with  the  origin  translated  along 
(T  -  t)w  for  each  t  in  the  interval  0  <  t  <  T,  and  v  -  v(t)  is  the  parachute  velocity 
vector  relative  to  air  In  the  horizontal  plana,  in  order  for  this  control  law  to  be  effective, 
the  initial  launch  position  must  be  in  the  vanity  of  the  free  descent  path  of  the  parachute 
which  would  hit  die  target  at  tcuchdov/n  In  the  presence  of  a  constant  known  wind. 
This  is  a  reasonable  assumption  in  view  of  the  limited  degree  of  control  available  in  the 
gliding  parachute  system. 

The  above  control  law  can  be  easily  updated  as  naw  estimates  of  the  wind  become 
available.  However,  it  Is  the  D.C.  or  mean  value  of  the  wind,  rather  than  the  instantaneous 
value,  which  should  be  used  as  the  estimate  since  the  control  law  is  optimal  only  for 
a  constant  wind  vector  w  on  the  ensuing  interval  T  >  t  >  to-  Although  (36)  is  only 
a  first  order  approximation,  the  term  In  brackets,  y,v2  -y2v,,  appears  to  be  a 
homogeneous  factor  in  each  of  the  terms  of  the  Taylor's  series  expansion  (27).  This 
was  only  validated  for  the  first  three  terms  of  (27)  and  applies  to  the  special  optimization 
problem  of  minimizing  the  performance  index  (5)  with  hi  ■  0  and  q2  >0. 

The  reciprocal  of  the  weighting  parameter  q2  in  (39)  represents  the  gain  of  the 
feedback  controller  which  must  be  adjusted  to  achieve  a  small  terminal  distance  error 
at  touchdown  without  incurring  excessive  bank  angles  during  the  parachute  maneuvering 
over  tha  interval  0  <  t  <  T.  The  banx  angle  <t>  upon  commands  from  the  controller 
electronics  must  be  designed  to  achieve  *ho  desired  relationship 

^  ■  tan-1  lViv2  -  y,v,]  (40) 

9  <h 

without  significant  delays. 

Future  work  should  include  the  following:  (i)  A  computer  simulation  of  the 
suboptima  I  control  law  and  a  comparison  of  the  results  with  optimal  control  trajectories 
obtained  by  solving  the  two  point  boundary  value  problem  for  various  initial  conditions. 
(I!)  Derivation  of  the  appropriate  filtering  equations  for  estimating  the  wind  based  on 
available  measurements.  Computer  simulation  of  the  controller-filter  combination  for 
random  wind  profiles.  (Hi)  A  formulation  of  the  problem  with  a  random  wind  which 
seeks  to  minimize  the  expected  value  of  the  terminal  distance  from  the  target,  in  this 
vein,  random  terminal  times  might  be  appropriately  studied  to  account  for  inaccurate 
altitude  data,  (iv)  A  formulation  of  the  problem  with  the  magnitude  of  the  parachute 
velocity  vector  a  control  variable  in  addition  to  the  directior. 
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Appendix  I 


Reduction  of  (18)  to  Elliptic  Integrals 

Successive  divisions  of  the  equations  in  (18)  by  x3'  yields  (with  €  defined  by  (21): 


dxj  _  cosx3  dx2  _  sinx3 

dx3  €x4  dx3  •  6x4 

dx4  x2(1)cosx3  -  xt(1)sinx3 

dx3  €  x4 


(41) 


The  third  of  these  equations  ccn  be  integrated  directly: 

xl  ■  C,  +  g  Ixjd)  cos  x3  -  J<a ( 1 )  sin  x3J  (42) 

where  the  constant  of  integration  C,  is  given  by 


Cm  -  [x4 ( 1 )) a — |r  *x3 ( 1 )  cos  x3(1)  -  x2(1)  sin  x3(1)) 


Using  the  relation  /  ^  *  -€/dr'  *  -€rfor  r»  1,  the  first  integral  is  obtained  as 

x4 

fv/c,  +  -  IFIk,*,)  -  F(k,  0O  )l  (44) 


where 

d  m\J |Xi(1)|a  *  |x2(1)|a  . 

ka«  — — — 

C,€+2d 

0O«  1/2  lx,(0)  +  tarr‘  JliHi  ) 

x,(1) 

0,-  1/2  [x3(1)  +  tan’1  1 

x,(1) 


(45) 


(46) 


and  F(M)  is  an  Elliptic  integral  of  the  first  kind. 


Substituting  (42)  into  the  first  of  the  equations  in  (41 )  leads  to  the  expression 
Xtd)  -  x, (0)  *  ± 


+  cos(x3(0)  +  0) 

•] 


cos(x3(1)  + 


i  &  wva  v 

^rrk 


(47) 


(1  —  ~ )(F(k, 0j)  -  F(k,0o) )  +  ~  (E(*,V»)  -  E(0o\aj)l 


where  4>0,  <Pi,  d,  C,  and  k  are  as  defined  above, 

0  -  tan"1  ,  sin*  a  <*  k* 

x,(1) 

and  E(0  Vs)  is  an  Elliptic  integral  of  the  second  kind. 

Substituting  (42)  into  the  second  equation  of  (41)  will  lead  to  a  similar  relation 
as  (47)  involving  Elliptic  integrals  of  the  first  and  second  kind.  This  equation  in  conjunction 
with  (44)  and  (47)  comprises  three  algebraic  equations  in  the  three  unknowns  x,  (1 ),  x3  (1 ) 
and  X4 (0),  given  the  initial  conditions  Xi(0),  x3(0)  and  x3(0j.  However,  the  arguments 
of  the  Elliptic  integrals  depend  on  the  boundary  data  in  a  mixed  way. 


Vector  Relations  for  Parachute  Dynamics  in  the  Horizontal  Plane 


Fig.  1 
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